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Abstract 

We generalize the pull approach to define the x 2 function to the analysis of the data with corre- 
lated statistical errors. We apply this method to the analysis of the Sudbury Neutrino Collaboration 
data obtained in the salt-phase. In the global analysis of all the solar neutrino and KamLAND 
data we find the best fit (minimum x 2 ) values of neutrino parameters to be tan 2 9\2 ~ 0.42 and 
<5m 2 2 ~ 7.1 x 10~ 5 eV 2 . We confirm that the maximal mixing is strongly disfavored while the 
bounds on <5m 2 2 are significantly strengthened. 
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As the solar neutrino physics moves from the discovery stage to the precision measure- 
ments stage increasingly more data are becoming available for a critical analysis. Several 
experimental groups recently announced new results. The Sudbury Neutrino Observatory 
(SNO) Collaboration announced results from the neutral-current measurements in the salt 
phase of the experiment A new, very precise measurement of the astrophysical S-factor 
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previous result [3|. 
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Earlier SNO charged-current (CC) measurements j^, |5[ confirmed the deficiency in the 
solar neutrino flux and helped narrow the neutrino parameter space. Earlier SNO neutral- 
current (NC) measurements [6( yielded a total (all flavors) 8 B solar neutrino flux which is 
in very good agreement with the Standard Solar Model (SSM) predictions 0. In a parallel 
development the reactor antineutrino disappearance observed by the KamLAND Collabora- 
tion [8( found to favor the same region of the neutrino parameter space (the so-called Large 
Mixing Angle (LMA) region) as the solar neutrino experiments. In a sense KamLAND and 
SNO are complementary experiments with KamLAND being more sensitive to the difference 
between squares of the neutrino masses, 5m 2 , and SNO providing increasingly more precise 
constraints on the mixing angle as well as 5m 2 . 

The new SNO NC measurement is able to determine the total active solar neutrino flux in 
a model-independent way. In particular no assumptions need to be made about the energy 
dependence of the flux. This feature makes the new data very valuable in restricting SSM 
or neutrino properties. 

In a previous paper P] we presented a global analysis of all the solar neutrino and Kam- 
LAND data using a covariance approach. The aim of this paper is to perform a similar 
analysis using the pull approach. We use the pull approach in the form proposed and im- 
plemented by Fogli et al Although the covariance and pull approaches are strictly 
equivalent, the pull approach seems to provide a treatment of statistical and systematic 
uncertainties which is more transparent, easier to implement, and computationally simpler. 
In our analysis we include correlations between statistical errors. 

We now enumerate the changes we incorporated in this paper as compared with our 
previous global analysis. The calculations of neutrino propagation in matter are done as 
described in Ref. j9|. Since the SNO Collaboration has not yet announced a live-time 
distribution for zenith angles, we used the live-time information previously provided in the 
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SNO data page 

[11] 

to calculate Earth regeneration effects and the 24-hour average neutrino 
survival probabilities. In our calculations we ignored experimental correlations between 
SNO's data sets as recommended by the SNO Collaboration They emphasized that 

even though treatment of these correlations require extensive knowledge of the SNO detector, 
they have very little impact on determining confidence level regions. 

The data set used in this analysis includes the average rate of the gallium experiments 
SAGE Q, GALLEX [14], GNO Q), the total rate of the Homestake chlorine experiment 
16 1, 44 data points from the 1496 days SK zenith-angle-spectrum n|, 34 data points from 



the 2002 SNO day-night-spectrum p as well as the newly-released SNO NC, CC, and 
electron scattering (ES) data 111. We calculated the CC cross section using the effective 



field approach 



18], with radiative corrections as described in Ref. I19J. In our calculations 



we used the counter-term value of L\& = 4.5 fm 3 



2l|. For new salt phase SNO data 



we calculated the detector-averaged cross-sections using updated detector response function 
given in Ref. Q] and higher threshold of 5.5 MeV as described in the SNO data page jl^ . 
The NC, CC, and ES fluxes are added as new data points. The systematic correlations of 
these new measured NC, CC, and ES fluxes are taken into account as described in [3]. The 
main systematic sources of errors, like energy scale, radial accuracy, isotropy mean, radial 
energy bias, internal background neutrons, neutron capture, etc. are taken into account 
with the appropriate correlation coefficients given in 12 1. 

Following the suggestion of Ref. [2^ , we added cross section errors first linearly and then 
quadratically in the parts contributing to the low and high energy parts of for radiochemical 
experiments. This change in cross section treatment resulted in a slight reduction of the 
allowed region contours. SNO collaboration has provided a statistical correlation matrix for 
the salt phase. The analysis of Ref. ^(| proves that for uncorrelated statistical errors the 
pull and covariance approaches yield identical \ 2 values. In the Appendix we write down 
the appropriate relations when statistical errors are correlated and prove that Fogli et al. 
result can be generalized to the correlated statistical errors. In our calculations we used the 
SSM systematics and fluxes Q]. 

We calculated the theoretical CC flux using the prescription given in the SNO data page 
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where 4>ssm is the SSM flux, P ee is the electron neutrino survival probability, and R is the 
energy response function. This result was used to calculate the CC/NC double ratio (see 



e.g. Ref. |23|). 

The collective effect of the fractional systematic errors of the SSM input parameters 
appear as a shift of the neutrino fluxes: 

*i - f&i, (2) 

where the index i represents the neutrino source (pp, 8 B, etc.). The shift fi can easily be 
calculated in the pull method (cf. Eq. (16) of Ref. 10]). In the covariance approach one 
can allow the 8 B flux to float freely. The shift in the pull approach provides a similar, but 
not identical function. 

All graphs are shown with several confidence levels. In the figures shown the darkest 
shaded areas are the 90 % confidence level regions. Lighter shaded areas progressively 
include the 95 %, 99 %, and 99.73 % confidence level regions. In this manner, for example 
the entire shaded region in a given graph is the 99.73 % confidence level region. Best fit 
points are marked by a dagger. Isolines are clearly marked with corresponding values of the 
quantities which are being examined. 

In order to compare the results before and after adding the SNO salt phase data, in 
Fig. [I]we show the LMA allowed region of solar neutrino parameter space. In this figure 
the CC/NC ratio isocountours, which are calculated for SNO experiment 2002 data with 
T e ,th = 5 MeV, are also shown. At best fit (marked by a cross), the value of this ratio is 
0.35. 

In Fig. 121 we show the allowed confidence levels in the neutrino parameter space using 
only the total NC, CC, and ES fluxes measured in the SNO salt phase measurement. This 
graph illustrates the crucial role of the additional spectrum information in constraining the 
neutrino parameters even after the neutral current measurements fix the 8 B flux. 

In Fig. El we show the allowed confidence levels in the LMA region of neutrino parameter 
space when chlorine, average rate of gallium experiments, SK zenith-spectrum and SNO 
salt phase measurement of NC, CC and ES data are used. In the calculations leading to 
this graph we excluded the SNO day-night spectrum data from 2002 in order to better 
understand the impact of new salt phase measurement. This figure can be compared to 
the Fig. [TJ Fig. 01 uses the same input information as Fig. [T] except the SNO day-night 
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information is replaced by the salt-phase results. Thus the relative impacts of the two data 
sets (salt-phase vs. day-night spectrum) in constraining the allowed neutrino parameter 
space are exhibited. Already the LOW solution is excluded at the 99.78 % confidence level. 
Clearly the information from the SNO salt phase helps significantly shrink the allowed region. 

The global analysis of all available solar neutrino data is illustrated in Fig. 0] (This 
figure is similar to Fig. El except that the SNO day- night spectrum data are also included). 
We find the best fit (minimum x 2 ) values of neutrino parameters to be tan 2 8u ~ 0.42 and 
5m\ 2 ~ 6.7 x 10~ 5 eV 2 from our combined analysis of all the solar neutrino data. Our 
minimum % 2 value is 70.5 for 83 data points and 2 parameters. The contours of the CC/NC 
ratio, which provides information about the electron neutrino survival probability, are also 
shown. Our best fit value for this ratio is 0.33. This is in good agreement with the SNO 
value of 0.306 ± 0.026(stat) ± 0.024(syst). 

The global analysis of all available solar neutrino plus the KamLAND data are shown in 
Figs. EJandlHl In Fig. Elthe 8 B flux is fixed at the SSM value whereas in Fig. |B]it is taken as 
a free parameter. In FigureElwe also show the ratio of this shifted 8 B flux to the SSM value. 
The best fit of the shifted flux corresponds to 1.02 times the SSM value. In both graphs 
we find the best fit (minimum % 2 ) values of neutrino parameters to be tan 2 # 12 ~ 0.42 and 
5m\ 2 ~ 7.1 x 10~ 5 eV 2 . Even though the allowed regions of the neutrino parameter space are 
quite similar in shape, allowing the 8 B flux float freely slightly enlarges the allowed region 
as one would expect. 

In conclusion one observes that the new SNO NC, CC, and ES data fluxes are statistically 
correlated since they are derived from a single fit to the data. In this paper we applied the 
pull method generalized to incorporate correlated statistical experimental errors to the anal- 
ysis of the SNO data. We confirm that the neutrino parameters are much better constrained 
after the addition of the SNO salt-phase data. In particular the maximal mixing is strongly 
disfavored while the bounds on 5m\ 2 are significantly strengthened. 
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Pull Approach with Correlated Statistical Errors 

Folllowing the notation of Fogli et al. we designate the statistical errors for the 
observable R n by u n , and the systematic error for the source k for this observable by c\. 
They also define scaled quantities 

A n = ^, (3) 

and 

q k n = t- (4) 
Thus for the covariance approach one can write 

Xcovar = E A n[pnm + E fc <M] _1 ^m, (5) 

n,m 

where p nm is the statistical covariance matrix including the correlated errors. For the pull 
approach one has 



Xpuii = min 

tSfe) 



(6) 



where is the univariate Gaussian random variable introduced in the pull approach. Min- 
imizing Eq. © with respect to we obtain 

ik = E ^ E PnL^mQn^ ( 7 ) 
h n, m 

where 

Sfch = [5kh + T,n, m < inPnm ( lmY 1 • ( 8 ) 

Defining the matrices S and D through their matrix elements 

B nm = J2lnC (9) 

and 

D nm = J2S hh qfot (10) 
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one can show that the covariance matrix is p + B. Using the definition given in Eq. (fTU|) it 
follows that the inverse of the covariance matrix is given by 

(p + B)- l = p- l -p- l Dp-\ (11) 

Substituting Eqs. (JJJ) and (fTTj) into Eqs. (0) and © after some algebra one gets 

Xcovar Xpull - (1^) 
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Ga CI SK SNO-Day-Night Data 
(CC/NC Double Ratio Isolines for SNO-Day-Night) 




tan 6 



FIG. 1: Allowed confidence levels in the LMA region of the neutrino parameter space when all 
solar neutrino experiments (chlorine, average of gallium, SNO and SK experiments) are included 
except the SNO salt phase data. CC/NC double ratio isolines for SNO experiment 2002 data also 
shown with T e t /j = 5 MeV. At best fit (marked by a cross), the value of this ratio is 0.35. 
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FIG. 2: Allowed confidence levels in the neutrino parameter space from SNO salt phase measure- 
ment of NC, CC and ES only. 
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SK Ga CI and SNO Salt Phase Data 
(SNO Day-Night Data Excluded) 
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FIG. 3: Allowed confidence levels in the LMA region of neutrino parameter space when chlorine, 
average rate of gallium experiments, SK zenith-spectrum and SNO salt phase measurement of NC, 
CC and ES included in the analyses. SNO day-night spectrum data from 2002 is not included 
in order to better understand the impact of the new salt phase measurement. Note that LOW 
solution is not allowed at 99.73% C.L. in this ^analysis (not shown). The best fit is marked by a 
cross. 



Ga CI SK SNO-Day-Night and SNO Salt Phase Data 
(CC/NC Double Ratio Isolines for SNO Salt Phase) 




FIG. 4: Allowed confidence levels of neutrino parameter space when all available solar neutrino 
data (chlorine, average gallium, SNO and SK spectrums and SNO salt phase data) included in the 
analyses. CC/NC double ratio isolines for SNO Salt Phase data are also shown with T et h = 5.5 
MeV. At best fit (marked by a cross), the value of this ratio is 0.33. 
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Ga CI SK SNO-Day-Night and SNO Salt Phase + KamLAND 
(Isolines for Ratio of Shifted 8 B Flux to SSM Value) 
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FIG. 5: Allowed confidence levels from the joint analysis of all available solar neutrino data 
(chlorine, average gallium, SNO and SK spectra and SNO salt phase) and KamLAND reactor 
data The isolines for ratio of the shifted 8 B flux (as described in the text) to the SSM value also 
superimposed on the plot. At best fit (marked by a cross) the value of this ratio is 1.02. 
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Ga CI SK SNO-Day-Night and SNO Salt Phase + KamLAND 

( 8 B Flux Free) 
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FIG. 6: Same as the previous figure, except 8 B flux is treated as a free parameter. Note that, 
shape of higher confidence levels look more similar to the original analyses of SNO Collaboration 
in which S B is treated as a free parameter as well. 
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